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ABSTRACT 

Sunyaev-Zel'dovich Effect (SZE) cluster surveys are anticipated to yield tight 
constraints on cosmological parameters such as the equation of state of dark en- 
ergy. In this paper, we study the impact of relativistic corrections of the thermal 
SZE on the cluster number counts expected from a cosmological model and thus, 
assuming that other cosmological parameters are known to high accuracies, on 
the determination of the w parameter and cxg from a SZE cluster survey, where 
w = p/ p with p the pressure and p the density of dark energy, and is the rms 
of the extrapolated linear density fluctuation smoothed over 8 Mpc/i~^. For the 
purpose of illustrating the effects of relativistic corrections, our analyses mainly 
focus on V = 353 GHz and Sum = 30 mJy, where v and 5*;^ are the observing 
frequency and the flux limit of a survey, respectively. These observing param- 
eters are relevant to the Planck survey. It is found that from two measurable 
quantities, the total number of SZE clusters and the number of clusters with 
redshift z > 0.5, and w can be determined to a level of ±1% and ±8%, re- 
spectively, with la uncertainties from a survey of 10000 deg^ Relativistic effects 
are important in determining the central values of and w. If we choose the 
two quantities calculated relativistically from the fiat cosmological model with 
(Js = 0.8284 and w = —0.75 as input, the derived as and w would be 0.819 and 
—0.81, respectively, if relativistic effects are wrongly neglected. The location of 
the resulting ag and w in the as — w plane is outside the 3a region around the 
real central cxg and w. 

Subject headings: cosmology: theory — galaxy: cluster — large-scale structure of 
universe 
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1. Introduction 

With fast advances in astronomical observations, the determination of cosmological pa- 
rameters has become one of the most important studies in cosmology (e.g., Eisenstein, Hu 
& Tegmark 1999; Spergel et al. 2003). Different kinds of observables are sensitive to dif- 
ferent underlying physical processes, and thus can be used to probe different cosmological 
parameters. The local abundance and the evolution of clusters of galaxies have been used 
extensively in the and erg determination, where Jl^ is the present matter density param- 
eter of the universe and ag is the rms of the extrapolated linear density fluctuation smoothed 
over 8 Mpc/i^^ (e.g., Bahcall & Fan 1998; Pen 1998; Haiman, Mohr & Holder 2001; Holder, 
Haiman & Mohr 2001; Fan & Chiueh 2001; Rosati, Borgani & Norman 2002). Supernova 
type la (SNela) (Schmidt et al. 1998; Perlmutter et al. 1999) and Cosmic Microwave Back- 
ground (CMB) radiation observations (e.g., Spergel et al. 2003) provide convincing evidences 
that we are living in a flat universe with about 70% of the matter being in the form of dark 
energy. As the properties of dark energy affect the global geometry of the universe and the 
growth of matter density fluctuations, the abundance and the evolution of clusters are de- 
pendent of the w parameter of the equation of state of dark energy, where w — p/p with p the 
pressure and p the density of dark energy (e.g., Caldwell, Dave, & Steinhardt 1998). Studies 
show that deep and large cluster surveys could constrain w to about 5% level (e.g., Mohr 
et al. 2002). In order to reach such a precision, however, we have to understand different 
types of systematic effects. Investigations indicate that the incompleteness of the knowledge 
about the intracluster gas (ICM) affects the cosmological parameter determination from X- 
ray cluster surveys, but influences less on the results from Sunyaev-Zel'dovich Effect (SZE) 
cluster surveys (Majumdar & Mohr 2003). In this paper, we study the impact of relativistic 
corrections of thermal SZE on the w and erg estimation from SZE cluster surveys. 

The typical gas temperature in clusters of galaxies is about a few kev, but clusters with 
temperatures as high as about 17 kev have been observed (e.g., Pointecouteau et al. 1999; 
2001). Thus relativistic corrections on SZE can be signiflcant (e.g., Birkinshaw 1999; Carl- 
strom. Holder & Reese 2002; Rephaeh 2002). For a SZE cluster survey, the selection function 
is mostly related to a mass threshold that depends on the sensitivity of the instrument and 
is also a function of redshift. For a survey with a relatively high flux limit, only clusters 
with large masses, and therefore high temperatures, can be detected. In this case, we ex- 
pect that relativistic corrections will impact the theoretical estimations of the survey results 
considerably. The amphtude of the change of SZE due to relativistic effects also depends on 
observing frequencies, with a relatively large correction around 400 GHz. For Planck, one 
of the best frequencies for SZE observation is 353 GHz and the flux limit is about 30 mJy 
(e.g., Diego et al. 2002). Both numbers indicate that we need to take relativistic effects 
into account in relating cosmological models to survey results. Diego, Hansen & Silk (2003) 
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discussed how the corrections affect the extraction of SZE signals from Planck survey. Here 
we study the impact on the SZE cluster number counts for different w models. Keeping 
other cosmological parameters fixed, we particularly investigate the influence on the w and 
CTg determination. 

For the redshift distribution of SZE clusters, the differences of different w models mainly 
show up at high redshifts. To extract w and erg, we then propose to use the total number of 
SZE clusters Nfot and the number of clusters with redshift z >= (denoted as Nz>z^) as 
constraints, where z^ is flux dependent. Thus to apply our method, one needs to know the 
redshift range of a cluster, but not its exact value. Therefore it provides an economic way 
to estimate cosmological parameters from large surveys, such as Planck. 

To demonstrate relativistic effects clearly, we consider single-frequency measurements 
with u = 353 GHz and the flux limit Sum = 30 mJy. Real Planck observations will be 
multi-frequency. Therefore for a particular cluster-detecting method, realistic simulations 
including survey characteristics are needed to quantify carefully the influence of relativistic 
corrections. More discussions on this regard are presented in Section 4. 

Our study flnds that for a survey of 10000 deg^, (Xg can be constrained to a narrow range 
with about ±1% accuracy around its central value if la uncertainties are allowed for both 
Ntot and iV2>o.5. The corresponding range for w is ±8%. Relativistic effects are signiflcant 
in terms of locating the central values of cosmological parameters. If they are ignored in 
deriving as and w from the observed Ntot and A^2>o.5, the misplaced location of {as,w) in 
as — w plane is outside the 3a region around the real central ones. 

The paper is organized as follows. Section 2 contains the relevant analytical formulation. 
Section 3 presents numerical results. Discussions are made in Section 4. 

2. Formulation 

2.1. Thermal SZE with Relativistic Corrections Included 

The thermal SZE, i.e., the thermal distortion of the CMB spectrum, is generated through 
the scattering of CMB photons with thermahzed hot electrons in, e.g., clusters of galaxies 
(e.g., Sunyaev & Zel'dovich 1980; Birkinshaw 1999; Carlstrom, Holder & Reese 2002). Be- 
cause the temperature of electrons is generally much higher than that of the CMB photons, 
inverse Compton scatterings are dominant, and the distribution of photons shifts toward 
high-frequency side. In terms of the dimcnsionless quantity x = hpV /ksTcMB^ where hp and 
ks are the Planck's constant and the Boltzmann's constant, respectively, and u and Tcmb 
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are, respectively, the frequency and the temperature of CMB photons, in the nonrelativistic 
hmit, the change of the CMB intensity can be written as 



AI 

— = Q{x)y, 

-'0 



where /q = {2hp/c^){kBTcMBf /h% 

Q{x) = 

and the Compton y-parameter 
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where rie is the number density of hot electrons, ctt is the Thomson cross section, Tgas is 
the intracluster gas temperature, rUe is the electron mass, and c is the speed of light. The 
integral is along the line of sight. 

In this limit, the shape of the distortion is independent of the gas temperature. For 
V < 218 GHz, Q{x) < 0, and the SZE shows as an absorber (or decrement) of CMB. For 
u > 218 GHz, Q{x) > 0, and the SZE shows as a source (increment) of radiation. 

For relativistic corrections, we use the analytical formula of Itoh, Kohyama & Nozawa 
(1998), which includes terms up to {kBTgas/rrieC^)^ . Then we have 

^ - 71^(^0 + 9eY, + 9lY, + 9% + 6%), (2) 

where 9^ = kBTgas/{meC^) and 

X 

° tanh(a;/2) ' 

and Yi, Y2, Y3 and Y4 are complicated functions of x [see Itoh et al. (1998) for detailed 
expressions] . 

In Figure 1 we show the single-scattering AI/Iq for ksTgas = 5 kev and for ksTgas = 
20 kev. For each pair, the solid and the dashed lines arc the results without and with 
relativistic corrections, respectively. It is seen that the effects arc significant for ksTgas = 
20 kev, with about 25% reduction of SZE signals at z/ = 353 GHz. It is also noted that the u 
value for the zero thermal SZE point depends on the gas temperature, in contrast with that 
of the nonrelativistic case (e.g., Birkinshaw 1999; Rephaeli 1995). As pointed out by Itoh et 
al. (1998), their perturbative approximation is valid up to ksTgas — 15 kev in comparison 
with that from exact numerical calculations. At ksTgas ~ 20 kev, the approximated result 
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deviates from the full numerical one significantly at u > 500 GHz. Atv < 353 GHz, however, 
equation (2) is still a good approximation even for ksTgas ~ 20 kev. 



Consider unresolved clusters, the total flux of the thermal SZE from a cluster at redshift 

z is 

^ AldA, (3) 



where Rd{z) is the angular diameter distance to the cluster, and the integral is over the 
surface area of the cluster. 

With the isothermal assumption, equation (3) can be written as 

^ = ^0 (ef_^ 1)2 (^0 + e,Y^ + elY2 + e% + etv,) y, (4) 



where 

Y 



If we denote as the gas fraction that is assumed to be a constant and X as the 
hydrogen mass fraction, then 

where M is the total mass of the cluster, and nip is the proton mass. Further under the 
equilibrium condition, the temperature Tgas and the mass M are related through the following 
relation (e.g., Wang & Steinhardt 1998) 

r 20^(7^ 1-3/2/ A/f N 2/3 

*Bi;„^0.944/,(l + .,|a„(0)AJ"3x[l-^] . (7) 

where — ff^/i/P with a constant factor of order unity that reflects the deviation from 
the simple spherical model, /i the mean molecular weight of the gas, and /3 the ratio of the 
kinetic energy of galaxies to the thermal energy of the gas. Qm{z) and ^a{z) are the matter 
density parameter and the dark energy density parameter at redshift z, respectively, and Ac 
is the ratio of the mean density of the cluster to the background matter density at redshift 
z. Here fi is related to X by /i = 4/ (5X + 3). 

Then given a flux limit for a SZE survey, one can get the corresponding mass limit 
Miim{Siim, z) from equations (4), (6) and (7). Note that the mass limit depends on cosmo- 
logical models. 
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2.2. Cluster Number Densities and Cosmological Models 



For the comoving number density of clusters, we use the fitting formula of Jenkins et 
al. (2001). Then we have 



^.,M) = 0.315^^1^1 exp[-|0.61 -ln(Z^.aM)|='-' 



dM 



Mau dM 



(8) 



where po is the present matter density of the universe, ctm is the rms of the linearly- 
extrapolated-to-present matter density fluctuation over the mass scale M, and is the 
linear growth factor of density perturbations. 



For (7m, we adopt the form given by Viana & Liddle (1999) 



R{M) 
8/i-iMpc 



(9) 



where h is the Hubble constant in units of 100 kms ^Mpc ^, the spatial scale R is related 
to M through 



R{M) = 



M 



and 



7(/?) = (0.3r + 0.2) 



(47r/3)po' 
2.92 + logio 



R 



.8/i-iMpc/J' 

with r the shape parameter of the linear density fluctuation spectrum. 
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In this paper, we consider flat cosmological models with the dark energy part described 
by the equation of state w — p/ p. The density parameter of dark energy is then specifically 
denoted by VLq. The linear growth factor Dz is given by the analytical formula of Wang & 
Steinhardt (1998) 

Dz^aexp{J^ ^[l-n^(zr]y (11) 



where a = 1/(1 + 2;) is the scale factor of the universe, and 

3 3 {l-w){l-3w/2) 



a 



5-w/{l-w) 125 (l-6w/5)3 



For the over-density of a cluster collapsed and virialized at redshift Zc, we have (Wang 
& Steinhardt 1998) 

3^1+5*a^3 



(12) 
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where Rta and Rmr are the radii of the cluster at the time of turn-around {z — Zta) and at 
the time of viriahzation {z — Zc), respectively, and 

Pcluster{Zta) /^TTX^ o.79+0.26Q„x-0.06m)| /i q\ 

^ " P,{zta) ~ ^''^ 

The redshifts 2;^^ and are related through t[zc) = 2t{zta), i.e., the time at Zc is twice the 
turn- around time, and 

Rta ~ 2 + r]t-3vv/2' ^ ' 

where 

_ ^^-i ^QjZta) 

Vt — Q / \ ) 

^^m\ Zta ) 

and 

-1/ + ZcY^q{zc) 



\ 1 2;t„ / 



3. Results of Analyses 

In the following studies, we take = 0.1, = 1.4, and X — 0.76 for the gas properties. 
For cosmological models, the Hubble constant Hq is taken to be 70 kms~^Mpc~^, Jl,n(0) = 
0.3, and O,q{0) — 0.7. For perturbations, we take F = 0.2. The observational frequency is 
taken tohe u — 353 GHz. 



3.1. Mass limits for Flux-limited SZE Surveys 

From equations (4), (6), and (7), it is clear that the corresponding mass hmit for a flux- 
limited SZE survey is cosmology-dependent, and this dependency is one of the factors that 
contribute to the differences of the expected SZE survey results for different cosmological 
models. 

In Figure 2, we show Mum, in units of IO^^H'^Mq, vs. redshift z for different w cos- 
mologies given Slim — 30 mJy. The three models are for w — —1, w — —0.75, and w — —0.3, 
respectively. For each pair, the lower one is the result without relativistic corrections, and 
the upper one is the one with the effects included. It is seen that the mass limit decreases 
as w increases. At 2; = 1, Miim{w = —0.3) ~ 0.7Miim{w = —1). The relative increase of the 
mass limit due to relativistic corrections ranges from 6.5% for w = — 1 to 7.5% for w = —0.3 
at z = 1. 
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3.2. SZE Cluster Number Counts 

As flat models with different w have different volume elements, linear growth factors for 
density perturbations, and mass limits, the redshift distribution of SZE clusters is sensitive 
to If, and thus can be used to probe the w value. 

In Figure 3a and Figure 3c, we show dN/{dQdz) for Sum = 5 mJy and Sum = 30 mJy, 
respectively, without relativistic corrections. In each flgure, results of three models with 
w = —1, —0.75 and —0.3 are plotted. For all the models, erg = 0.85. Therefore the figures 
show the effects of different w on the cluster redshift distribution. We see that given a 
flux limit, models with larger w (smaller absolute value of w) predict more SZE clusters 
over almost the whole redshift range considered. In Figure 3b and Figure 3d, relativistic 
corrections on dN/ (dQdz) are shown. It is seen that they are larger at higher redshifts, and 
relatively smaller for higher w at redshifts z >1. For Sum = 30 mJy, the relative changes are 
about 13% at z = 0.5, and at 2; = 2, they are about 41%, 38% and 32% for w — —1, —0.75 
and —0.3, respectively. 

In Figure 4, we show the surface number density dN/dfl vs. w. For each pair, the upper 
and the lower ones are for the results without and with relativistic corrections, respectively. 
For Slim = 30 mJy, the respective dN/dQ = 1.45 and 1.30 at w = —0.75, with the change 
of about 10%. For Planck, the survey area is about 10000 deg^. Then the corresponding la 
uncertainty of dN/dfl due to Poisson fluctuations is about 0.012. Thus the change due to 
relativistic effects is about 12a. This indicates the significance of the effects in extracting 
w value from observations. For example, with single-frequency measurements, given the 
observed surface number density of 1.30, one would infer w ~ 0.9 if relativistic corrections 
are neglected, which is about 20% off the real value oiw — —0.75. We elaborate this further 
in the next subsection. 

In Figure 3 and Figure 4, erg = 0.85 for all the models. That is, we assume that erg, 
as well as other cosmological parameters, have been well determined from other observa- 
tional information. On the other hand, studies have shown that the total number of SZE 
clusters is sensitive to ag, thus potentially, as and w can be constrained from SZE sur- 
veys simultaneously. From the shape of the redshift distributions in Figure 3, we see that 
with the same total number of SZE clusters, the differences of the models are still appar- 
ent. For Slim = 5 mJy, the number of clusters with 2; > 1 is significantly larger for larger 
w. The ratio dN{z > l,w = -0.3)/ d^ to dN{z > l,w = -l)/dn is about 1.8. For 
Slim — 30 mJy, the number of clusters with z > 0.5 shows strong dependence on w, with 
[dN{z > 0.5, w = -0.3) /dn]/[dN{z > 0.5, w = -l)/dQ] ^ 1.8. These analyses suggest that 
it is possible to constrain both ug and w using the total number of SZE clusters and the 
number of SZE clusters with z > z^, with z^, flux-limit dependent. With this methodology. 
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we only need to know the redshift range of a cluster, but not its precise z. Thus it is rela- 
tively easy to be realized in practice. In the next subsection, we show how well we can apply 
this method to constrain erg and w. We also discuss relativistic effects on the parameter 
determination. 



3.3. Constraints on erg and w 

We concentrate on extracting information on erg and w from SZE cluster surveys by 
assuming that other cosmological parameters have been prc-dctcrmined. The values of those 
other parameters considered in this paper were listed at the beginning of Section 3. 

Two pieces of information will be used to constrain erg and w: the total number of SZE 
clusters and the number of clusters with z > Zm- Here we take Sum — 30 mJy and Zm — 0.5. 

The model with erg = 0.8284 and w — —0.75 is adopted as the fiducial one. We choose 
w — —0.75 (but not w — —1) as our central value so that relativistic effects can be seen 
clearly (see Figure 6 in the following). Our analyzing procedures are as follows. Given the 
total surface number density dN/dfl calculated from the fiducial model, we search for erg 
for different models such that they have the same dN/dfl as that of the given value. In 
this way, we get a relation between erg and w. Similarly, there is another crg-w relation from 
the number density of clusters with z > 0.5 [denoted as dN{z > 0.5) /dQ]. The two lines 
from the two relations should intersect at erg = 0.8284 and w = —0.75. Given a survey area, 
we can estimate the possible constrained ranges for erg and w allowing different uncertainty 
levels. 

Figure 5 shows ag-w relations for Sum = 30 mJy without (thin lines) and with (thick 
lines) relativistic corrections. The set of relatively flat lines are from dN/dQ, and the other 
set of fines are from dN{z > 0.5) /dfl. The solid lines are from the fiducial dN/dO, and 
dN{z > 0.5) /dfl, respectively. The dashed fines are ±3er results for dN/dfl, and the dotted 
lines are ±3er results for dN{z > 0.5) /dfl. Here Poisson statistics is applied to estimate er, 
and the survey area is taken to be 10000 deg^. 

For the non-relativistic case, we find that the sofid fines can well be fitted by 

as{w) = 0.8284|w + 0.75 - i|(o.i353-o.oi50|«,+o.75-i|)^ ^^5^ 

and 

as{w) = 0.8284|«; + 0.75 - i|(o.i809-o.oo56h+o.75-i|)_ (^g) 

To avoid crowding, we did not plot ler and 2er lines in Figure 5. Our studies show that with 
ler uncertainties in both the total number of SZE clusters and the number of clusters with 



-10- 



z > 0.5, (78 and w are constrained to be within the ranges (0.8209, 0.8360) and (—0.81, —0.69), 
respectively. Thus for a survey such as Planck, the respective accuracies of the determined 
(Tg and w can possibly reach to a level of about ±1% and ±8%. The 3a determinations are 
about ±3% for erg and about ±25% for w if other cosmological parameters have already been 
known to good precisions. 

In the relativistic case, the approximated analytical relations between erg and w for the 
solid lines are 

a,{w) = ag(-0.75)|w + 0.75 - i|(o.i246-o.oii5i«,+o.75-i|)_ ^^^^ 

and 

as{w) = as{-0.75)\w + 0.75 - i|(o.i6969-o.oooo9|«,+o.75-i|)_ ^^g) 

The relations are shghtly flatter than those of the nonrelativistic case (eq. [15] and eq. [16]). 
As seen in Figure 3d, relativistic effects are stronger at higher redshifts. It is also known 
that models with larger w (lower absolute w value) predict more high redshift clusters, and 
thus the relativistic corrections to clN/dQ and dN{z > 0.5) /dQ are larger. This results the 
decrease of the differences between different w models in terms of the two measurements, 
and therefore the flatter o"g-w relations. The la and 3a ranges for cxg are (0.8213,0.8356) 
and (0.8056,0.8491), respectively, and the respective ranges for w are (-0.815,-0.69) and 
(—0.95, —0.57). We see that relativistic corrections do not affect the precisions of the pa- 
rameter determination very much. 

On the other hand, however, if relativistic effects were neglected in extracting ug and 
w from a SZE survey, the determined central values can be significantly off the real ones. 
To estimate this misplacement, we take the relativistic dN/dO, and dN{z > 0.5) /dQ from 
the fiducial model with erg = 0.8284 and w — —0.75 as inputs, but search for as-w relations 
with nonrelativistic mass limits. In Figure 6, we plot the results as solid lines. The contours 
are the la (solid), 2a (dotted) and 3a (dashed) lines around cxg = 0.8284 and w = —0.75 
(denoted by the star). Note that the solid contour is composed by ±la lines for the total 
number of SZE clusters and the ±la lines for the number of clusters with z > 0.5. The 2a 
and 3a contours are similar. First let us consider one parameter case (given the other) using 
the total number of SZE clusters as the constraint. Given cxg = 0.8284, the derived central 
w from the wrong solid line is w ^ 0.92, while the 3a value is w 0.78. On the other hand, 
iiw — —0.75 is given, the derived as is as ~ 0.814, while the 3a value is as ~ 0.824. We see 
that the misplacement is dramatically large. If both erg and w are considered as unknown, 
the derived values from the cross-point of the two solid lines are o"g ^ 0.819 and w ~ —0.81. 
Thus the derivations from the flducial values are about 1.1% for cg and 8% for w. Although 
this offset is smaller than those of the one parameter cases, the derived "central values" are 
outside the 3a region (the region surrounded by the dashed contour). Therefore relativistic 
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effects are very significant and cannot be ignored. 

4. Discussion 

As tlie most important property of SZE is its redsliift-independence, SZE cluster sur- 
veys are one of tfie best to probe tfie structure formation at high redshifts. Because of the 
cosmological model dependence of the geometry of the universe, of the evolution of den- 
sity fluctuations, and of the mass limit corresponding to a given flux limit, the redshift 
distribution of SZE clusters is sensitive to cosmological parameters. 

We particularly discussed the effects of relativistic corrections of thermal SZE on the 
redshift distribution of clusters, and further on as and w determinations from SZE cluster 
surveys. The corrections affect the cluster redshift distribution through changing the mass 
limit for a given flux limit. The changes depend on the observing frequency and on the flux 
limit. For u = 353 GHz and Sum = 30 mJy, at 2; = 1, the relative increase of the mass limit 
due to relativistic effects is about 6.5% and 7.5% for w = —1 and w = —0.3, respectively. 
With as = 0.85 for all the models, the corresponding decrease of dN/dVt is about 8.9% and 
13.5% for w = —I and w = —0.3, respectively. For dN{z > 0.5) /dfl, the percentages of the 
decrease are, respectively, about 18.4% and 21.3% for w = —1 and w = —0.3. 

With fixed other cosmological parameters, we studied the constraints on as and w from 
the total number of SZE clusters and from the number of clusters with z > Zm, where Zm 
is flux-limit dependent. For Sum ~ 5 mJy, Zm ~ 1-0 is appropriate. For Sum — 30 mJy, 
Zjn — 0.5. The nice part of our methodology discussed here is that we do not need to 
know the precise redshift of a cluster, but only its range, i.e., larger or smaller than z^- 
Therefore it is applicable in large surveys, such as Planck, because then one can use not- 
so-precise methods to estimate the redshifts of clusters, such as the photometric method 
or morphological redshift estimates (Diego ct al. 2003). Our study showed that at Sum = 
30 mJy, for a survey of 10000 deg^, cg can be determined to a level of ±1% and w to an 
accuracy of ±8% with la uncertainties in both the total number of clusters and the number 
of high redshift clusters {z > 0.5 in our analyses). With 3a uncertainties, the determined 
parameter ranges are ±3% and ±25% around the central values of ag and w, respectively. 
In our uncertainty estimations in subsection 3.3, we assumed that the redshift information is 
available for 100% SZE clusters. That is, when we assessed the a for the number of clusters 
with z > 0.5, we used the full number of clusters with z > 0.5 for a survey of 10000 degl If 
the redshifts are available for 10% of clusters, with la uncertainties, the determined ag and 
w have the ranges (0.815, 0.842) and (—0.89, —0.62) around the central values as = 0.8284 
and w = —0.75. Thus with 10% samples having redshift information, the accuracies of la 
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determination for as and w are ±1.6% and ±18%, respectively. Relativistic corrections do 
not have significant effects in tfiis aspect. 

However, we found tfiat the central values of as and w can be significantly misplaced if 
relativistic corrections are not included in extracting cosmological information from observed 
cluster number counts. For Sum — 30 mJy, the offset for as is about 1.1%, and about 8% for 
w. In terms of the percentages, this misplacement is comparable to the la deviation. But 
the location of the wrong central values in the a^ — w plane is significantly outside the 3a 
region around the real central {as,w) as seen in Figure 6. Therefore, relativistic effects are 
important, and should be considered carefully in analyzing survey results. 

We need to emphasize that our discussions on relativistic effects concentrate on single- 
frequency measurements at 353 GHz. In reality, Planck will conduct multi-frequency obser- 
vations, and the overall relativistic effects could be weaker than the case of u = 353 GHz 
depending on the techniques used to detect SZE clusters. Diego, Hansen & Silk (2003) stud- 
ied the impact of relativistic corrections on the reconstruction of the y map from simulated 
AT/T of different frequency channels with Bayesian non-parametric method. They found 
that for clusters of T = 10 kev, the relative difference in the recovered y between relativistic 
and non-relativistic approaches in the component separation was about 4% for low-redshift 
clusters and about 8% for high-redshift clusters. This difference is smaller than that at 
z/ = 353 GHz, which is about 15%. This is because the y maps were constructed through the 
weighted average of AT/T of different frequency channels, and relativistic corrections are 
smaller than 10% for u < 353 GHz. Thus the single-frequency estimation of relativistic ef- 
fects presented in this paper can be an over-estimation. To quantify the impact realistically, 
one can, with simulated maps including characteristics of a survey, get two sets of mass lim- 
its of detectable SZE clusters at different redshifts with and without relativistic corrections 
for a particular cosmological model. Then analytical studies similar to those shown in this 
paper can be carried out based on the two sets of mass limits. We would like to mention 
that the fiux limit for completeness, thus the mass limit, depends on the method applied to 
find clusters. The limit Sum = 30 mJy (at u = 353 GHz) adopted in our analyses is from 
the estimation of Diego et al. (2002) based on the Maximum Entropy Method (Hobson et 
al. 1998; Hobson et al. 1999). On the other hand, Bayesian non-parametric method gives 
the completeness limit ~ 200 mJy (at u — 353 GHz) (Diego et al. 2003), and a multi- 
filter approach of Herranz et al. (2002) finds Sum ~ 170 mJy. With such high fiux hmits, 
relativistic corrections can be large even with the multi-frequency diminishing effect taken 
into account. 

The constraints on as and w shown in the paper are derived under the assumption that 
all other cosmological parameters have been determined by other studies. If we allow, e.g.. 
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Jl^(O), to vary as well, the estimation can be affected considerably (e.g., Haiman, Holder, 
& Mohr 2001). On the other hand, however, there are indeed different types of observations 
that can give good constraints on those parameters. For example, galaxy surveys, such 
as two degree field Galaxy Redshift Survey (2dFGRS) (Colless et al. 2001) and the Sloan 
Digital Sky Survey (SDSS), can put constraints on the shape of the perturbation spectrum, 
which primarily depends on flmh (Szalay et al. 2001; Dodelson et al. 2002; Percival et al. 
2001; 2002). Combined with the measurements on the Hubble constant Hq, fi^ can be well 
constrained. In any case, studies on the multiple-parameter estimation with the methodology 
put forward in this paper and on the impact of relativistic corrections are desired, and will 
be carried in the future. 

We thank the referee for the helpful comments and suggestions. This research was 
supported in part by the National Science Foundation of China under grant 10243006, and 
by the Ministry of Science and Technology of China under grant TGI 999075401. 
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Fig. 1. — Thermal SZE AI / after a single scattering with Jo = {'2hp/c'^){kBTcMB)^/hp 
vs. frequency u for ksTgas = 5 kev and 20 kev. In each pair, the solid line represents the 
result without relativistic corrections, and the dashed line is for the results with relativistic 
corrections up to {kBTgas/meC^)^. 
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Fig. 2. — The mass limit in units of W^^h~^MQ vs. redshift z. The flux limit is 

Slim = 30 mJy. The solid lines are for w = —1, the dotted lines are for w = —0.75 and 
the dashed lines are for w = —0.3. In each pair, the upper one corresponds to the one with 
relativistic corrections included, and the lower one is for nonrelativistic M^^. 
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Fig. 3. — The redshift distribution of SZE clusters dN/ (dQdz) in units of deg~^ for — 
5 mJy. No relativistic corrections are included, and as — 0.S5. The solid line is for w = — 1, 
the dotted line is for w = —0.75, and the dashed hne is for w = —0.3. 
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Fig. 4. — Relative changes of the redshift distribution of SZE clusters due to relativistic 
corrections for Sum = 5 mJy and Ug = 0.85. 
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Fig. 5.- 



Same as fig. 3a but for Sum — 30 mJy 
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6. — Same as fig.3b but for Sum — 30 mJy 
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Fig. 7. — The surface number density of SZE cluster dN/dQ in units of deg~^ vs. w. 
The sohd, dotted, dashed and dash-dotted hues are for Sum = 5 mJy, 30 mJy, 60 mJy and 
90 mJy, respectively. In each pair, the upper one is the nonrelativistic result, and the lower 
one is the result with relativistic corrections included. Here cxg = 0.85 for all the models. 
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Fig. 8. — The as — w relations without (thin hnes) and with (thick hnes) relativistic cor- 
rections. The set of relatively flatter lines are from the total number of SZE clusters, and 
the steeper lines are from the number of clusters with z > 0.5. The solid lines are from the 
central values of Nfot and Nz>o.5, respectively. The dashed lines arc the ag — w relations from 
Ntot =t 3cr, respectively, and the dotted lines are the respective relations from Nz>o,5 ± 3a. 
The survey area is taken to be 10000 deg^. 
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Fig. 9. — The misplacement of the central and w. The two sohd hnes are from the 
nonrelativistic calculations given Nf^t and Nz>o,5 calculated relativistically with w — —0.75 
and as — 0.8284 (represented by the star in the figure). The survey area is taken to be 
10000 deg^. The sohd, dotted, and dashed closed hnes are la, 2a and 3a deviations around 
w — —0.75 and as — 0.8284 (see the text for details). 



